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In center vortex theory, beyond the simplest picture of confinement several conceptual problems 
arise that are the subject of this paper. Recall that confinement arises through configuration aver- 
aging of phase factors associated with the gauge center group, raised to powers depending on the 
total Gauss link number of a vortex ensemble with a given Wilson loop. The simplest approach 
to confinement counts this link number by counting the number of vortices, considered in d=3 as 
infinitely-long closed self-avoiding random walks of fixed step length, piercing any surface spanning 
the Wilson loop. Problems arise because a given vortex may pierce a given spanning surface several 
times without being linked or without contributing a non-trivial phase factor, or it may contribute a 
non-trivial phase factor appropriate to a smaller number of pierce points. We estimate the dilution 
factor a, due to these inert or partially-inert vortices, that reduces the ratio of fundamental string 
tension Kf to the areal density p of vortices from the ratio given by elementary approaches and 
find a = 0.6 ± 0.1. Then we show how inert vortices resolve the problem that the link number of a 
given vortex- Wilson loop configuration is the same for any spanning surface of whatever area, yet a 
unique area (of a minimal surface) appears in the area law. Third, we discuss semi-quantitatively 
a configuration of two distinct Wilson loops separated by a variable distance, and show how inert 
vortices govern the transition between two possible forms of the area law (one at small loop separa- 
tion, the other at large), and point out the different behaviors in SU(2) and higher groups, notably 
SU(3). The result is a finite-range Van der Waals force between the two loops. Finally, in a problem 
related to the double-loop problem, we argue that the analogs of inert vortices do not affect the fact 
that in the SU(3) baryonic area law, the mesonic string tension appears. 

PACS numbers: 11.15.-q, 12.38.-t, 11.15.Tk UCLA/06/TEP/01 



To some extent, our understanding of area laws in confining gauge theories is based on intuition and plausibility. 
Certainly, there can be no doubt that, in the fundamental representation of SU(N), the expectation value (W) of the 
trace of a simple flat Wilson loop T involves the area A of the flat surface spanning it and not the area of any other 
spanning surface: 



[Here Kp is the fundamental string tension, and we use imaginary anti-Hermitean gauge potentials, incorporating 
the gauge coupling g in them.] But, at least in the center vortex picture of confinement, it is not always easy to see 
how some of these plausible results follow from the basically simple mechanism of confinement, based on linkages of 
vortices with Wilson loops. In this paper we discuss several confinement puzzles, all of them connected by the theme 
of inert vortices. By this we mean vortices that do not link to Wilson loops in the usual way, but occupy space that 
could have been occupied by truly-linked vortices. For brevity we also use this term to refer to partially-inert vortices 
that arc linked, but with a smaller link number than would naively be expected. 

We seek the effects of inert vortices only for large Wilson loops, those whose length scales are all large compared 
to the fundamental gauge-theory length A. As the loop scales approach A the effects of inert vortices either disappear 
or are substantially modified. 

In the center vortex picture, the area law arises through group-center phase factors raised to powers depending 
on the Gauss link number of the vortex condensate with the Wilson loop. This link number can be calculated from 
the intersections of a vortex with any surface spanning the Wilson loop. To characterize the condensate we will 
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stick for simplicity to d=3, although there is no real qualitative difference between three and four dimensions. In 
d=3, vortices are closed stringlike tubes of chromomagnetic flux, a finite fraction of which have infinite length, and 
in d=4 they are closed 2-surfaces, whose description raises complications. In d = 3 we model the vortices as closed 
self-avoiding infinite-length random walks on a cubical lattice of lattice length A. This model is similar in spirit, if 
not in implementation, to the d = 4 models for SU(2) and SU(3) center vortices given by others [l|, for use in lattice 
computations. It is also similar to the usual identification of center vortices as infinitesimally-thick objects, called 
P-vortices 0, that live on a lattice dual to the lattice where Wilson loops live. 

Even for the archetypical example of a flat Wilson loop for SU(2) 7 inert vortices are important. In this simple case, 
there are many types of inert vortices, including those that pierce the surface twice within a few characteristic lengths 
A, or that pierce it an odd number of times but have link number less than the number of pierce points. The net 
effect of inert vortices for a simple flat Wilson loop is that the density of linkage differs from the density of piercing 
by a factor a, which we call the dilution factor, lying between zero and one. This factor is a rough but useful estimate 
of the various ways in which vortices can be inert. In Sec. IIIII we estimate that a lies in the range 0.6±0.1. In the 
dilute gas approximation (DGA), the usual result for the SU(2) fundamental string tension Kf is Kf — 2p; dilution 
by inert vortices modifies this to 

K F = 2ap. (2) 

[For SU(3), the standard DGA gives K F = 3p/2; the diluted DGA is K F = 3a 3 p/2, as for SU{2). It may be that the 
ST/ (3) dilution factor a 3 is not precisely the same numerically as it is for 5/7(2), but our estimates are not accurate 
enough to see much of a difference. For general SU(N) there are in principle as many dilution factors as vortex 
densities.] 

It would, of course, be good to compare our estimates with lattice data. Unfortunately, it has turned out 0,0,0 to 
be rather difficult to calculate the density p on the lattice, for a number of reasons. Among these are the dependence 
on gauge of the center-vortex location procedures; effects of Gribov copies; and finite-size effects. Ref. Q states that 
lattice artifacts are so important that these authors cannot really find a reliable value for p. However, p| claims a 
value of Kf / p of 1.4, which taken literally might indicate a dilution factor around 0.7, at least if the DGA is more 
or less correct. The best way to attack the numerical estimation of dilution might be to simulate directly a model of 
self-avoiding random walks, in the spirit of 0, rather than to work with QCD itself. 

A second major issue arises because link numbers can be calculated (through Stokes' theorem) by counting intersec- 
tions of vortices with a surface spanning the Wilson loop, but any spanning surface can be used for the link-number 
calculation, not just the surface that ultimately appears in the area law. So it is not clear what area is to be used in 
the area law, nor even why there is a unique area. We discuss these issues in Sec. IIVI showing how unlinked vortices 
resolve the paradox of two or more possible areas for a simple Wilson loop. 

The third issue, elaborated in Sec. |V] is an interesting variant on the question of when vortices are linked or not. 
This issue was raised before at a qualitative level |(J. We consider two identical Wilson loops separated by a certain 
distance, and ask how the overall VEV of these two loops depends on separation. For SU{2) this is a problem 
somewhat like the corresponding soap-bubble problem, where there are (at least) two minimal soap films that can 
appear for two wire frames close to one another. We can think of each loop as a qq meson, and our results indicate 
a Van der Waals potential between the two loops, which breaks (as also happens for soap bubbles) at a critical 
separation between the loops. Similar but more elaborate results hold for ST/(3). 

Finally, in Sec. I VII we consider the baryonic area law for St/(3). This has been explored in center vortex theory 0, 
where it is shown that the area law comes from three surfaces with quark world lines and a central line as boundaries. 
To some extent the baryonic area law seems to have issues like those of the double Wilson loop in Sec. [3 in particular, 
there might appear to be correlations like the Van der Waals potential of the two-loop problem that could modify the 
accepted baryonic area law. However, we show that analogs of the inert vortices in the two-loop problem, which are 
in fact not inert even though they are in the same geometry except for orientation effects, do not affect the fact that 
the linearly-rising potential for the three quarks has precisely the mesonic string tension for each of its three sheets. 
Sec. IVIll contains a summary. 

Except for Sees. ITU IlIII we use the DGA approximation to describe our results. But it turns out to be convenient, 
for these two other sections, to use a standard form of the area law that contains the DGA as a limiting low-density 
case. 

II. BASIC PRINCIPLES FOR THE CENTER- VORTEX AREA LAW 

The center- vortex picture for gauge group SU (N) invokes a vortex condensate. If N > 3 there are N — 1 types 
of vortices labeled by an integer fc, 1 < k < N — 1 that gives the vortex magnetic flux in units of 2ir/N. A vortex 
labeled k is the antivortex of the vortex labeled N —k. These vortices are characterized by an areal density pk — PN-k 
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for each vortex type. By this areal density we mean (in all dimensions, not just d = 3) that the average number of 
/c-vortices that pierce any flat surface of area A is p^A. There is not much theoretical insight into the values of these 
different densities for SU(N) with TV > 3. However, for SU(2) and SU(3) there is only one density, which we term 
p, that sets the scale for the string tension (in SU(3) there are two types of vortices, but one is the antivortex of the 
other and they have the same densities). In this paper we will only consider gauge groups SU(2) and SU(3). 

A condensate of vortices can form only if a finite fraction of them has essentially unbounded length (or an unbounded 
number of steps in the random walks describing the vortices). Only such vortices, long compared to any Wilson loop 
scale, can contribute to the area law. These can be linked or not, depending on the circumstances we encounter. 
Finite- length vortices are therefore inert, in our terminology, and will be accounted for by a renormalization factor 
that we will not attempt to calculate here. 

Consider now a simple flat Wilson loop and the flat surface spanning it. This surface is taken to lie in a plan of 
the lattice dual to the vortex lattice and is divided into squares of this dual lattice; we call these A-squares. Any such 
square is pierced by a single vortex with probability p. This probability is related to the areal density of the vortex 
condensate by 

P = P X 2 ; (3) 

the probability that a square is unoccupied is p = 1 — p. As on the lattice, p can be extracted 0, 0, 0, Ej from the 
VEV of a square Wilson loop one lattice unit on a side. Denote this VEV as (W(l x 1)); then for SU(2) 

(W(lxl))=p-p; p=i[l-(W(lxl))]. (4) 

The only difference from the lattice definition is that in lattice computations the length scale for the Wilson loop is 
not a physical quantity A but a lattice spacing, and the lattice version of p must be scaled via the renormalization 
group to find a physical probability, such as used in the present model [Eq. ©]. Note that p is bounded by 1/2, 
since (W(l x I)) lies between and 1. Moreover, note that the probability p as derived from Eq. Q is not subject 
to dilution, which applies only to large Wilson loops. 

The assumption that vortices in different A-squares are statistically independent leads to the standard argument 
for center-vortex confinement, which ignores inert vortices. The confining area law (discarding perimeter effects) for a 
Wilson loop follows from an ensemble average of center-group elements. Each of these is of the form Yii Zf ki , where 
Zi is an element of the center of the gauge group as specified by the properties of the i th vortex (and the group 
representation of the loop itself). For example, for the fundamental Wilson loop in SU (2) the only non-trivial element 
of the center has Zi = — 1. In the above product, Lki is the Gauss linking number of this vortex with the Wilson 
loop. The Gauss linking number, a topological quantity, can be written as an intersection number of the vortex and 
any surface spanning the Wilson loop; its (integral) value is independent of the choice of surface. In the SU(2) case 
the necessary average is 

(exp[in^2 Lki}) . (5) 

i 

Aside from the assumption of independent A-squares, the critical assumption for expressing confinement in the 
center- vortex picture is that p is the probability that a vortex is actually linked once to a flat Wilson loop. When an 
odd number of vortices is linked once, the Wilson loop has value -1 and when an even number is linked, the value 
is +1. If the assumption is true, the area law follows from multiplying the probabilities p — p of Eq. for all the 
A-squares of the spanning surface. 

Another useful way of expressing this area law is to write out the combinatorics for vortex occupancy of N$ sites 
of a surface spanning a Wilson loop T: 

(Wr> - P Ns - Nsp^p + Ns ( N s-V p Ns-y + ... = [p- p ]Ns = [i _ 2 p]^A 2 . (6) 

Here the number Ns 3> 1 of sites on a given spanning surface S is Ns = As/ A 2 , where the surface has area As, and 
each term represents the number of ways of arranging empty and once-filled A-squares. 

Next we need to modify the area law for dilution, which arises from several factors: If a vortex penetrates an even 
number 2N P of times, there are N p sites that lead to unit phase factor in the Wilson-loop VEV, although these sites 
are occupied. In effect, the vortices filling such sites are inert (although they may, strictly speaking, be linked to the 
Wilson loop topologically) . Similarly, if a vortex penetrates an odd number (greater than I) of times it is linked and 
gives a non-trivial phase factor, but three or more sites are occupied, rather than the single site assumed when we 
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related the string tension and the piercing probability as in Eq. ©. We will argue in Sec. IIIIBl that the diluted form 
of the standard equation 10 is 

(W r ) = {l-2ap} A s/ x \ (7) 

yielding a string tension 

K F = ~ln(l-2ap). (8) 

The DGA approximation is the small-p limit of either Eq. J^J (undiluted) or of Eq. JHJ) [diluted; sec Eq. Note that 
while the effect of dilution on the probability p is simply to renormalize it, the effect of dilution on a dimensionless 
quantity such as Kp/p cannot be characterized as a renormalization of a dimensionful quantity such as the density p 
or the A-square area A 2 : 

_L = ln(l-2a;)U=apA». (9) 

p x 

It seems very difficult to resolve dilution problems completely by analytic methods; the best we can do is to give a 
semi-quantitative discussion of how these factors renormalize downward the linkage probability from p. 



III. THE DILUTION FACTOR FOR A FLAT SURFACE 



In this section, dealing with the gauge group SU(2), we distinguish between the previously-introduced probability 
p that a vortex pierces a A-square, thus contributing unit link number, and the density of link number per A-square, 
which is what we really need. This density is reduced by inert-vortex effects. We attempt to capture, in some mean- 
field sense, these effects approximately by introducing a dilution factor a that effectively reduces the pierce probability 
p to ap. 

We begin by classifying various ways in which vortices can pierce a spanning surface yet not be linked (in the sense 
that they are associated with a trivial phase factor), or are linked but are to be associated with a reduced density of 
linkage numbers. For brevity we refer to all these as inert vortices. 



A. Types of inert vortices 



It is useful to distinguish three types of inert vortices; only Types II and III need detailed discussion. Type I vortices 
have finite length, and correspond in some sense to localized particles. The vortex condensate may have some of these, 
but they cannot explain confinement, since for large Wilson loops those that are linked contribute only to sub-area 
effects such as perimeter terms in the VEV. Presumably the effect of such vortices is essentially to renormalize the 
areal density of vortices of unbounded length. We will not discuss such finite- length vortices any further, so from now 
on we are only concerned with the issue of vortices much longer than any Wilson loop scale, and the extent to which 
these are or are not inert. 

The distinction between the remaining two types of long vortices is this. Type II vortices exhibit what we will 
call local return, by which we mean that a vortex, however long, penetrating a localized flat surface has its highest 
probability of returning to that surface after only a few more steps of the random walk. This probability is not to be 
confused with the probability that a random walk will be near to where it has been after a large number N s of steps; 
this probability decreases like N^ d ' 2 , where d is the dimension of space-time (see the Appendix). For the present 
section, dealing with a flat spanning surface, this type is important, but once the renormalization of p has been made, 
they are not very important in the further applications of Sees. II VI M 

Finally, a Type III vortex is one which, having penetrated a surface, penetrates it a second time with high probability 
after a large number of steps. This can only be true, in view of the above remarks, if the surface is of a special type, 
including those used in Sees. HVllVl and I VII such as a closed surface or one that has curvature radii comparable to its 
length scales. An infinite-length vortex must penetrate a closed surface at least twice. 

A qualitative model of Type III effects might be to assume that after renormalization for the effects of Types I and 
II vortices, the remaining effects are seen for vortices composed of infinitely-long straight lines intersecting surfaces 
compounded of flat segments. The rationale for the straight-line vortices is that the local returns have been accounted 
for by renormalization for Type II vortices. In this model no vortex can return locally to a single surface element and 
so explicit Type II vortices are missing; more than one flat surface segment must be involved, and these must form a 
non-planar surface. We will not use this model in the present paper. 



(a) 



(b) 




(c) 



(d) 



FIG. 1: In these figures, the surface being pierced lies in a plane on a lattice (shown) dual to the vortex lattice (not shown), 
(a) Simple piercing by a vortex (thick line), (b) An inert vortex, penetrating the surface with two more steps after the first 
piercing, (c) An inert vortex that has rendered the square labeled X inaccessible to another vortex, (d) A partially-inert vortex, 
using up three A-squares for a single linking. 

Among the Type II vortices, the only ones to be considered in this section, there is essentially an infinite number 
of subtypes. Several of them are shown in Fig. d 

The first, Fig. ^a), shows the standard penetration assumed in the usual area-law formula of Eq. ©: Unit link 
number associated with a single piercing. Fig. ^b) shows an inert vortex (zero link number) produced by two 
additional steps from the first piercing. Fig. dc) shows another inert vortex, using three additional steps, that 
makes the square labeled X inaccessible for piercing by another vortex, because of the mutual- avoidance requirement. 
Fig. \\^d) shows a vortex with four extra steps that is linked, but has pierced three A-squares, thereby again reducing 
the A-squares available for more vortices. Note that in every case shown in the figure except for Fig. ^a) the link 
number is reduced relative to the piercing number. 

There are other ways for vortices to be inert, for example, a vortex may have a link number that is a multiple of N 
for SU(N) with N > 2. Related effects take place in the baryonic area law for SU(3) (Sec. IVII . 



The fiat surface is divided, as before, into Ns A-squares. Of these, on the average pN$ are pierced once by a vortex, 
and of the pierced squares, a fraction a on the average contribute to the area-law formula with a minus sign (in 



We use a simple statistical model, ignoring shape effects and certain correlations and assuming that dilution is 
statistically independent of piercing. It is easy to see that the coefficient of 



B. 



Definition of a 



5(7(2)). 




(10) 



in the formal expansion of 



1 = \p + ap + (1 



a)p\ Ns 



(11) 
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is the statistical weight for a configuration with K2 + K3 pierced A-squares, of which K2 are going to give a minus 
sign in the Wilson loop VEV. We then have 

(W) = [p- ap + (1- a)p] Ns = [1 - 2ap] Ns , (12) 

showing that, as previously specified, a simply renormalizes p. For the DGA formula for the Wilson-loop VEV we 
find 

(Wr> = exp[-2apA 8 ], (13) 

in which the string tension is renormalized by the factor a from its previous value. The DGA is perhaps made more 
plausible by dilution. 

Not all the effects of inert vortices can be captured by the simple dilution factor defined above. Vortices are 
correlated with each other through self-avoidance, as in Fig. ^c), and the specific geometry of the portion of a 
random vortex walk that penetrates the spanning surface more than once can matter. Unfortunately, even simpler 
problems cannot be solved analytically. For example, the statistics of co-existing but mutually- and self-avoiding 
monomers and dimers (see [8| for a mean-field approach and earlier references), has no exact solution (except in the 
limit of close-packed dimers |2j|). This is because dimers and higher multimers do not obey simple (e. g., multinomial) 
statistics. One can appreciate this from the observation that on an empty lattice of N$ sites a dimer can be laid 
down in 2N$ ways, but only N$/2 self-avoiding dimers can be put down in total. (For monomers, of course, these 
two numbers are the same, namely N$.) 

In the inert-vortex problem there are, in principle, a huge number of multimers, of various shapes and sizes, that 
should be accounted for. Rather than attempting some elaborate generalization of the monomer-dimer problem, we 
proceed as follows. 



C. Estimating the dilution factor 



The problem is to estimate the probability that a self-avoiding random walk of infinite length, having pierced the 
Wilson surface (which we call W) once, pierces it again one or more times. We need two different types of probabilities 
for this problem. The first is the standard probability density p^(N; rh) that a random walk on the d = 3 lattice is at 
the lattice point rh after exactly N steps. If the difficult restriction of self-avoidance is dropped, this is given by (see, 
e. g., E3) 



l 3 r 2lT 

p 3 {N;m) = ,o 0/v [TT / <%][cos0i + cos6 2 + cos0 3 ] N exp[i6- m]. (14) 



This is normalized so that 



p 3 (N = 0;m) = 5 iA x (15) 



and the sum over all rh of this density yields unity. In the limit of large N and components of rh (which is an 
integer- valued vector) p 3 (N;rh) has the usual Gaussian form. This is not the probability of real interest, although 
it can be used in certain circumstances to find the probability that we need. We will call the one that we do need 
q 3 (N; W; rh), the probability density that a random walk piercing at the origin re-pierces for the first time at rh using 
exactly N steps, with the vector rh restricted to lie in the lattice plane that is nearest neighbor to the Wilson surface. 
The probability density P3(N; W; rh) is simply P3(N; rh) with rh restricted to the Wilson surface. Our interest is in 
probabilities summed over the Wilson surface, so we define 

p 3 (N;W)= P^N-W-rh); q 3 (N;W) = £ q 3 (N;W;rh). (16) 



By the standard rules of probability |ll| 



N 

p 3 (N; W) = 6 Nt0 + ]T <? 3 (J; W)p 3 (N - J- W) (17) 



which says that the random walk, having pierced the surface for the first time after J steps, may penetrate it many 
times again before ending on the surface after N steps. [Note that it takes N = 2 additional steps of the random 
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walk, at minimum, to re-pierce the surface, given that these steps are counted as starting with the first step after the 
first piercing.] If the total number of piercings is odd (even) the vortex is linked (unlinked). The case of no further 
intersections after the first piercing is to be included; the standard area-law model of Sec. [H] is equivalent to assuming 
that this probability is unity, and all other probabilities are zero. This is far from the case, as we will see. 
Eq. (|17|l is easily solved in terms of generating functions 

P 3 (s;W)=J2P3(N;W)s N ; Q 3 (s;W) = ^q a (N;W)8 N . (18) 

2 

We have 

P 3 (s;W) = 1 + P 3 (s;W)Q 3 (s;W); Q 3 (s;W) = l- / — . (19) 

P 3 (s;W) 

Note that we may express the elementary solution for P 3 (s; W), which is 

p ^ w) = T^kwy < 20) 

in a suggestive way for the original probabilities: 

P3 (N; W) = 1 + E W) + £ Q3(K; W)} 2 + ... (21) 

showing how the probability p 3 is compounded from probabilities of first return [see, e. g., Fig. [Jd)]. 

The final probability of interest is the probability that the random walk ever re-pierces the surface; this is clearly 
given by the sum of all the q 3 , or by Q 3 (s — 1; W). 

Equation i|19|l also holds in certain other problems, notably the gambler's ruin problem (described in the Appendix), 
which asks for the probability that a d = 1 random walk starting at the origin ever returns to it. In this problem 
backtracking (re-tracing the last step) is allowed, and (as we review in the Appendix) it is straightforward to show 
that the corresponding probability of ever returning is exactly unity. It is also unity in two dimensions. However, our 
problem differs from the gambler's ruin problem in two essential ways: No backtracking is allowed, and ultimately the 
problem becomes three-dimensional, when the radius of gyration iV 1 / 2 of the random walk becomes large compared 
to the size scale L of the Wilson loop. [Actually, the radius of gyration for self-avoiding walks has an exponent 
somewhat different from 1/2, but we ignore that complication here.] It has long been known (as reviewed in |10|) 
that the probability of return to the origin in d = 3 is finite, with a value of approximately 0.34. Furthermore, 
self-avoidance completely changes the problem; for example, a self-avoiding "random" walk in d = 1 has probability 
zero of ever returning to any site it has reached. 

Let V be the probability that a vortex, having penetrated the Wilson surface once, never penetrates it again. Then, 
with Q 3 (\\ W) as the probability that it ever penetrates the surface again, 

V=l-Q 3 {V,W). (22) 

Or one may write, as in Eq. (|21|l . a probability sum rule saying that the sum of probabilities of piercing exactly once, 
exactly twice, etc., must be unity: 

i = ni+E w) + E q ^ K; w ^ 2 + ■ ■ ■>■ ( 23 ) 

As before, the sum over K begins with 2 for self-avoiding random walks. 

It is now necessary to estimate V, or equivalently Q 3 (1;W). Unfortunately, this is not a straightforward matter 
when there is self-avoidance (see the Appendix for a brief review of some of the well-known analysis when this condition 
is not imposed). One way to proceed is simply to count the number of self-avoiding paths going from an original 
piercing to another piercing as a function of their step length, and look for ways of partially re-summing the results. 
Our proposal, given below, is more accurate for random walks that do not backtrack than for true self-avoiding walks; 
the difference is that a non-backtracking walk may violate the condition of self-avoidance by looping back on itself. 
We propose the following expression for non-backtracking walks, valid for J > 2: 

B(J!H0 _ ( * ) 1 [( | + 1,« + (|_»^. (24, 

The corresponding generating function, approximately valid for self-avoiding walks, is 
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The explanation of the terms is as follows. The factor 4/25 is the probability for Fig. QJb), for J — 2. For a non- 
backtracking walk in d — 3 there are 5 possible choices to add a new step to the random walk. The horizontal step 
in this figure has probability 4/5 (note that we are summing over all possible sites), and the next, vertical, step has 
probability 1/5. The appearance at larger J of 3/5 in the formula expresses the probability that a non-backtracking 
random walk will take its next step horizontally (with respect to the Wilson surface), and 1/5 is the probability of a 
vertical step in one particular direction (up or down). [The restriction to even powers of 1/5 is easily understood by 
drawing a few figures, in the style of Fig. ^] So a random walk once started in a horizontal plane has a three times 
larger probability of staying in that plane than of moving to a plane higher above the Wilson surface. It is easy to 
check the combinatorics of the low-order powers of 1/5 in Eq. (|24() . and also the highest powers, for example, the last 
term (4/25)(l/5) J ~ 2 . This is the probability that the random walk goes as high as possible, takes one horizontal step, 
and then returns straight down to the Wilson surface. This can only happen for even J, and one easily sees that the 
maximum attainable height is J/2. We have compared the approximation of Eq. (|24|l with explicit counting of self- 
avoiding walks through J = 8 and find that the difference between self-avoiding and non-backtracking is acceptably 
small. 

The difference between non-backtracking walks and self- avoiding walks first appears at N — 5, where there are 
144 non-backtracking walks but only 128 self-avoiding walks. This is about an 11% error, but because the erroneous 
contribution to the N — 5 walks is 16/3125 ~ 0.0052, out of a total of about 0.533, the error in the final result from 
N = 5 terms is less than 1%. In any case, to deal with true self-avoidance rather than just non-backtracking raises 
the same issues as having a condensate of mutually-avoiding vortices, and we do not discuss that issue in any detail 
in this paper. 

We find an approximation to Qs(l; W) by summing over all J in Eq. (|24|l . with the result 

Qsd; W) * (±) x I[ T -1 M + T -L M ] = 1 * 0.533. (26) 

At this point one might guess that the dilution factor a should be equal to the probability V of never returning . 
But this is not quite right; we really need the probability of being linked or unlinked. To do this, one should separately 
find the probabilities that there is an even or odd number of piercings, as expressed by separating the even and odd 
powers of Qz{\) in Eq. l)2.'iH . We use Vl for the probability that a vortex is linked, and Vjj for the probability that it 
is unlinked; these are defined by 

even t \ / t \ / 

Now it seems that a = Vl- If we to use our approximation <5s(l) = 0.53 we would find the link probability (dilution 
factor) to be about 0.65, somewhat larger because we are counting three, five, . . . piercings as well as a single piercing. 

Eq. I|27fl is not quite right either, because while Vl expresses the probability of linkage, when (say) three sites are 
used for the link instead of one, there is dilution that must be accounted for. We account for this by calculating a sum 
of weighted probabilities, where the weight for the appearance of K powers of Qz(\) in Eq. I|27|) is the inverse of the 
number of sites occupied by the linked vortex, or (K + 1) . For example, when K = 2 three A-squares are pierced, 
as shown in Fig.^d). This result for a, when modified by the original piercing probability p, should give something 
like the link number per unit area, in effect increasing the size of A-squares to account dilution. The weighted sum 
for Vl gives another estimate for a: 

^ N + 1 * 20 3 1 L l-Q a (lV 



even 



For Q 3 (l) = 0.53 this gives a = 0.52. 

It is plausible that Qz(\;W) lies between the no-self-avoidance d — 3 value of 0.34 and the d = 1 value of 1. In 
that case, the estimate of a with no inverse-site weighting [Eq. (|27j gives a finite value even in the d = 1 limit, where 
we find Vl — Vjj = 1/2, so the dilution factor is 1/2. This is a singular limit, because the original probability V 
used to construct Vli'Pu vanishes, but this is cancelled by a singularity in the sum over Qs. It is not surprising 
that the limiting probabilities are each 1/2, since there is no way of distinguishing an even number of piercings from 
an odd number. The value of a from this equation for the d = 3 value Qz{t]W) = 0.34 is about 0.75. Eq. J2HJ>, 
with inverse-site weighting, gives 0.75 for Qs(l;W) = 0.34 and for Q3(1;W) = 1. Presumably this latter case is 
unrealistic, because at some point the radius of gyration of the random walk is large compared to the Wilson loop 
scale L and the problem really is three-dimensional; in any case, the d = 1 self-avoiding case is completely opposite 
to the no-self-avoidance case. 



FIG. 2: A flat Wilson loop, and two spanning surfaces, labeled £ and S. 



To be more accurate in estimating a, one would need to account for mutual avoidance effects such as shown in 
Fie. [Tl c) . and be more precise about weighting various random- walk configurations. We will not attempt that here, 
and close by saying that our estimates for the dilution factor are consistent with a — 0.6 ± .1. 

In the following sections we look for further manifestations of inert vortices of Type III, going beyond the local 
effects that led to the dilution factor. But these local effects still occur, and so everywhere in our arguments the 
original probability p should be replaced by the diluted probabilities 

p = ap; p = 1 — p. (29) 



IV. INERT VORTICES AND THE SPANNING SURFACE 



Now we return to the problem that the spanning surface used for counting link numbers as piercings is arbitrary, 
yet there surely is a unique area in the area law. The simplest possible case, where there is no doubt as to the answer, 
is that of a flat Wilson loop, and one can get the right area law by using the flat spanning surface in Stokes' theorem. 
Yet one should also get the right answer by using any other surface. How does this come about? 

Fig. (J2J shows such a flat Wilson loop, with two spanning surfaces. The first, labeled E, is flat and is correct. The 
second is labeled S. We choose orientations so that the combined surface E + S is oriented. Of course, this surface is 
closed. 

We must now improve upon the techniques outlined in Sec. IIVI to account for inert vortices of type III. The 
calculation of the area law based on the flat surface E needs no change. But what if we instead wished to calculate 
the area law based on surface S 1 ? Here there is extra dilution. The vortices linked to surface S are still those linked to 
surface E, which must pass through surface E and surface S once each, but there are also inert vortices, which pass 
through surface S twice and surface E not at all. [We will not account for vortices linking three, five, . . . times, most 
of which are included in the dilution factor.] We denote the number of linked vortices by Nl and the number of inert 
vortices by Nj/2. The factor of 1/2 in the latter definition simply reflects the fact that every inert vortex pierces 
surface S twice, so that Nj is the total number of pierce points of inert vortices. It is not convenient to introduce this 
factor of 1/2 for linked vortices, because the diluted pierce probability p, introduced above, is related to the number 
Nl of vortices linked to surface E by 

N L = ^ (30) 

where is the area of surface E. That is, Nl is the number of pierce points of vortices on surface E. 
The total number of vortex piercings of the combined surface E + S, denoted N to t, is 



N tot = 2N L + Ni, 



(31) 



10 



with the two arising because linked vortices penetrate both surface S and surface E. By hypothesis of a uniform areal 
density of vortices, this total number of vortex pierce points on the combined surface E + S is also given by 



Combining these equations, one finds 



N tot = i^±^l. (32) 



= Kjg-ZM . (33) 
A 



To calculate the area law in the DGA, we can easily use the formulas of Sec. [H] for the flat E surface, which has a 
total of JV S = A^/X 2 A-squares: 

(W) = [p - p] Ws — e - KpAs , (34) 

with Kp = 2p/X 2 as before. It is more interesting to calculate it from the point of view of the other spanning surface 
5. Here we must account for the diminished probability of occupation of this surface by linked vortices, since some of 
them, as counted by Njj, are inert. So we change the probabilities by adding to p, the probability of no occupation, 
the probability pi = Nj/Ns of occupation by an inert vortex of Type III. This gives 

p->p =p + pi = p+p{l -—) = 1 — . (35) 

As As 



Similarly, the new link probability p = 1 — pQ is 



p=— — . (36) 
As 



Using Eq. (|34|) with the new probabilities for the surface 5 one gets, going to the DGA limit as before: 

(W) = [po - P] Ns -> e- 2 ^ As = e - W 2 Ms . (37) 



V. COMPOUND WILSON LOOPS 



Consider the compound Wilson loop shown W comp in Fig. [21 composed of two identical but oppositely-oriented 
Wilson loops. [The orientation is irrelevant in SU(2).] They are separated transversely by a distance z. The relevant 
expectation value is 

(W comp ) = (W(l)Wt(2)) (38) 

When the distance z between the loops is small compared to the loop, intuition |6j suggests that the surface whose 
area should appear in the area law, as in Eq. is the minimal surface 3 joining the two loops as shown in Fig.|3jA). 
When z is large, intuition suggests that the situation in Fig.|3{B) holds, where each loop is spanned by its own minimal 
surface with no connection to the other surface. We will show that intuition is indeed correct for the compound area 
law of center vortex theory, and give an approximate interpolation formula for intermediate values of z. 

Begin with the SU{2) case, where the orientation of the Wilson loop does not matter, and W = . As usual, we 
assume that the time extent T, the spatial extent R of the Wilson loop, and their separation z are all large compared 
to the QCD scale length A and assume that T R. In this limit the spanning surface S3 of Fig.|3jA) is nearly flat, 
and has area A3 ~ 2zT. (We ignore the contribution ~ 2zR from the top and bottom.) The Wilson loop surfaces 
5*1,2 have areas A\ = A2 = RT. 

Imagine now a configuration where all surfaces 5*1.2,3 exist, so that there is a closed surface with two marked 
contours, the Wilson loops 1 and 2. This constitutes a minor generalization of the configuration already considered in 
Sec . 11 VI There are several ways that vortices can be linked or inert (in the sense of Sec. lIII A"|l . after the renormalization 
of Type I and Type II vortices. Use the notation iVj for the total number of vortices penetrating surface Si. These 
obey 



(39) 
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(A) (B) 

FIG. 3: Two identical but oppositely-oriented rectangular Wilson loops 1 and 2, separated by a distance z. (A) When z is 
small compared to loop scales, the Wilson-loop area, labeled 3, connects the loops. (B) When z is large, the Wilson-loop areas 
are the disjoint areas, labeled 1 and 2, of each loop. 



Each of these is subdivided as follows: The number of vortices piercing S% and S3 an odd number of times is called 
jVi3, with analogous notation for A^3 = Ni 3 ; the number piercing surface Si and S2 each an odd number of times is 
N12] and the number entering S3 and then exiting the same surface is -/V33. 

As indicated in Sec. the probability that a vortex known to be at a point remote from a surface that then 
penetrates the surface (once) is proportional to the surface area, and is finite when the QCD scale length A vanishes. 
If we require the number of vortices penetrating one surface, say Si, and then another surface, say S2, that probability 
is proportional to A\A 2 . So we now assume that 

N 12 ~A 1 A 2 = Al; N w ~AiAx N 33 ~ A\. (40) 

These vortex-linking numbers are related to the total vortex-piercing numbers by 

Nx = N u + Nxa\ N 3 = 2N 13 + N 33 . (41) 

Eqs. (|39I40I41|1 are easily solved to yield 

2pA\ . _ pAiA 3 . pA\ 
Nl2 ~ \*(2A 1+ A 3 y Nl3 ~ \i(2A 1+ A 3 y ~ X 2 (2Ai + A3)' (42) 

We will now compute the expectation value (W) of the product W = W1W2 of the two Wilson loops, from the 
point of view of the surfaces Si 2- Only the number N 13 contributes non-trivially to an SU(2) Wilson loop. As in 
Sec. HVI we introduce modified probabilities 

. _ A 2 Ar 13 pA 3 . 

Just as in calculating the Wilson loop VEV from Eq. <|37|l . we have in the DGA: 

(W comp ) = (ft, -p)W - exp-[^^] = ex P -[^^]. (44) 



This formula is only approximate, but it shows features that we believe are generally correct (and one feature 
that is not correct). For example, the heavy-quark potential V [coefficient of — T in the exponent of Eq. 144|l ] has 
the behavior V ~ 2Kpz in the limit z <C R, showing that the the area law of two SU (2) Wilson loops [or, as in 
Fig. |21 two oppositely-directed Wilson loops] disappears as the two loops approach each other and form an AT-ality 
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FIG. 4: A baryonic Wilson loop in 5(7(3) is composed of three simple Wilson loops sharing a common central line (expanded 
in the figure). The central line is invisible to SU(3) center vortices. 



zero configuration. (For iV-ality zero loops there is a pseudo-area law, coming from the finite size of vortices [Ilj , at 
distances z ~ A, but that is irrelevant here.) In the opposite limit of z ^> R we find V — 2KpR, so the VEV is just the 
product of the separate VEVs for two Wilson loops. While this in itself is correct, the approach to this limit cannot 
be, for it would yield a residual potential V — 2KpR which vanishes only like 1/z. In actuality, at some point when 
R ~ z the spanning surface switches on the length scale A from S3 to Si + S2, much as the corresponding soap-bubble 
surface would for two physical wire loops 1 and 2. We do not see this breaking because we have not included effects 
coming from a network of gluon world lines running from loop 1 to loop 1, from loop 2 to loop 2, and from loop 1 
to loop 2. The simplest step in the formation of this network has been discussed [l3| in connection with baryonic 
and mesonic hybrids, having extra gluons along with their quark content. There it is shown that a single gluon acts 
like a physical string which separates a Wilson-loop spanning surface into two surfaces, and the string requires extra 
energy to be stretched if the string stretching leaves the minimal spanning surface. The fluctuations associated with 
this stretching should yield the Luscher [3] term, as well as leading to the transition from surface S3 to Si + S2 as 
the gluons in the network recombine. But when this stretching is not too great the potential should behave roughly 
as given in Eq. 144|) . 

In SU(N) with N > 2 the orientation of the Wilson loops does matter. If the loops have opposite orientation, as 
shown in Fig. [31 the problem is essentially the same as for SU(2). But one should also consider the problem of two 
loops of the same orientation. Since for SU(3) the treatment of this problem is quite similar to that for the baryon, 
which is a compound of three Wilson loops, we defer further discussion to Sec. IVII 



VI. VORTICES ANALOGOUS TO INERT VORTICES AND THE BARYONIC AREA LAW 

It is by now well-established both in center vortex theory and on the lattice 0,^1 that, as shown in Fig.0] in 
SU(3) the heavy-quark baryonic potential has three Wilson-loop surfaces spanning a boundary consisting of the quark 
world lines, with the three lines coinciding along a central line (the so-called Y law). This is a three- fold compound 
Wilson loop, sharing some of the features of the double loop discussed in Sec. [V] In particular, a vortex may pierce 
two loops. Such a vortex is not inert, as were the vortices piercing surfaces 1 and 2 of Fig. |5]and discussed in the 
previous section. Instead, linkage with these analogs of inert vortices are important to establish that the string tension 
in the baryonic area law is precisely the mesonic string tension Kp. 

The analysis of the dilution factor a for SU (3) is slightly more complicated, because a vortex piercing twice is 
not unlinked as it is in SU(2); it is simply equivalent to an antivortex piercing once. Given the complete symmetry 
between vortex and antivortex, this means that an antivortex piercing twice is equivalent to a vortex piercing once. 
There is a modification of the numerical value of the dilution factor, since it is now possible to dilute unit link number 
by occupying a minimum of two sites, rather than a minimum of three as for SU(2). However, the string tension is 
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FIG. 5: We introduce a new surface labeled 4 that creates a closed surface consisting of the two Wilson-loop surfaces for 
quarks 1 and 3 plus a surface joining these. 



still lessened by a single dilution factor a, just as for £{7(2), and we can continue to use the notation developed for 
that case. We will not pursue the question of what the value of the dilution factor is for SU(3); it must be quite 
similar to that for 5(7(2). 

The SU(3) version of the standard area law for a single Wilson loop, given in Eq. © for SU(2), is 

(W) = {p+ |[e 2 ™/ 3 + e -2W3]}A s /A 2 = [1 _ M]A S /A 2 _ (45) 

Here we use p/2 as the probability (accounting for dilution) that a vortex of flux 27r/3 pierces the surface, with equal 
probability that the antivortex of flux — 27r/3 pierces it; this means that, as before, p = 1 — p. This leads immediately 
to the DGA SU(3) string tension K F = 3p\- 2 /2. 

To discuss the baryonic area law in similar terms to those used with the double Wilson loop, introduce a surface 
made of flat pieces, labeled 4 and shown in Fig. that forms a closed surface when combined with the surfaces 
spanning Wilson loops 1 and 3. 

In what follows we will be careful to distinguish vortices from antivortices, although at the end the symmetry 
between them makes this simply a distinction for convenience of exposition. As in Sec. we introduce the number 
JVi3 as the number of times that the same vortex (not antivortex) pierces both spanning surface 1 and surface 3 an 
odd number of times, and AT 14 is the number of times a vortex pierces both surface 1 and 4. We consider that the 
dilution factor a has been applied, so we can take this odd number to be just one. The number Ni is the number of 
vortex piercings on surface 1, and it obeys 

JVi = ^ = N w + N u . (46) 

The factor of 2 in the denominator arises because there are just as many antivortices piercing any given surface; it is 
the same factor of two dividing the explicit probabilities in the area law of Eq. (|45|l . 

Any vortex that pierces surface 1 once must leave through either surface 3 or surface 4. If it leaves through surface 
4 it is associated with a phase factor exp[27ri/3]; if it leaves through surface 3 it is associated with a phase factor 
exp[— 2m/3]. For any such vortex configuration there is an equally probable antivortex configuration, with the phase 
factors interchanged. So the total number of piercing, vortices plus antivortices, is 2N\, and the number of antivortices 
piercing surface 3 is N13, etc. This, plus the relation of Eq. I|4bl . means that the standard calculation of the baryonic 
area law using surface 1 (plus surfaces 2 and 3), giving the Y law, is exactly the same as accounting for the linkages 
of vortices and antivortices through surfaces 1 and 4, plus linkages through surfaces 1 and 3 (plus permutations for 
the other two quarks). There really are no unlinked Type III vortices in this problem, the way there were in Sees. II Vl 
Ivl The baryonic string tension is precisely the mesonic string tension, because when it is calculated from the three 
quark surfaces 1, 2, and 3 separately, without regard to the other two quarks, it is just the standard simple Wilson 
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loop calculation. But it is just the same when linkages to two quarks are considered. Of course, no one doubted the 
equality of string tensions for baryons and mesons, but the point was to see how it worked out in the center vortex 
picture. 

VII. SUMMARY AND CONCLUSIONS 

We have discussed semi-quantitatively some of the effects of inert vortices that do not couple as effectively as they 
might to a Wilson loop, and that change the probability of linkage by a dilution factor a < 1. In turn, this changes the 
SU(2) DGA string tension from Kp — 2p to Kp — 2ap, where p is the areal density of vortices piercing a large flat 
surface. Our estimate, based on partial sums of non-backtracking walks, is a — 0.6 ±0.1. Understanding inert vortices 
also leads to understanding of how it is that there is a unique area in the area law, even though any surface spanning 
a Wilson loop is suitable for counting link numbers in the center vortex picture. We found a finite-range Van der 
Waals force, due to inert vortices, that links two mesonic Wilson loops, in this case vortices that are simultaneously 
linked to both Wilson loops. And finally we showed that the analogs of inert vortices in the two-loop problem do not 
interfere with the usual formulation of the area law for SU (3) baryons, based on considering the three Wilson loops 
as independent. 

There seems to be no possibility of a detailed analytic approach to these problems, which therefore are best studied 
further with lattice computations. One can, of course, create a center vortex condensate through simulation of the 
underlying non-Abelian gauge theory, but it would also be very interesting and possibly simpler to study inert-vortex 
effects with simulations that begin with an a priori vortex model, similar in spirit to the approach of Qj, rather than 
having to deal with all the complications of the full gauge theory. Indeed, it seems that relating Kp and p on the 
lattice has numerous complications 0. [3- 0] ■ and has not been attempted recently. But the simulation of two Wilson 
loops seems approachable on the lattice. 
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APPENDIX A: PROBABILITIES ISSUES 



This appendix discusses two issues: 1) A brief review of well-known analysis for return probabilities for random 
walks with backtracking allowed; 2) The probability for a vortex to pierce two separated Wilson loops. The first serves 
for cautionary notes in the problem calculating the dilution factor of Sec. lIIII and the second arises in finding the area 
law for the compound Wilson loop of Sec. [V] For the most part and for dimension d < 4, analytic results are only 
available for random walks without self- avoidance constraints, and we will only discuss those in the next subsection. 



1. The gambler's ruin problem [Sec. ITTT| 

The only problems that one has any real chance to analyze are for non-self- avoiding walks, and we will review a 
famous one, the gambler's ruin problem, here. The basic probability concepts are the same for self-avoiding walks, 
and so this review should be helpful to those unfamiliar with the underlying ideas. The cautionary note here is that 
while the d = 1 gambler's ruin problem seems superficially quite similar to those of Sec. IIIII in practice they are very 
different, both because the questions of Sec. lIIII have three-dimensional effects and because they deal with self-avoiding 
walks. 

We ask for the probability qi(K; m) that an unbiased random walk starting from the origin on a lattice of points 
in d = 1 will return to the point m > for the first time after K steps. The case m = is the standard gambler's 
ruin problem. For an unbiased walk (all step probabilities equal to l/(2d) in d dimensions) Polya long ago proved 
that in the limit of infinite steps the probability of return to the origin is unity in d = 1,2 but less than one in all 
other dimensions. This is perhaps plausible from the fact that the probability of being at the origin after N steps in 
d dimensions is, for N ^> 1, 



Pd 



(N) ~ N- d / 2 , (Al) 



and the sum over N diverges at large N in d = 1, 2 but not in higher dimensions. 
As in Sec. IIIII the relation between the p±(N; m) and the qi(N; m) is 

N 

Pl (N;m) = J2Pi( N - J;m = 0)q 1 (J;m) { Pl (N = 0;m) = S mfi ). (A2) 
j=i 

Note that in the sum on the right p\ has m set equal to zero, because we are compounding the probability of first 
return to m, as given by q\, with the probability of being at this same point m after more steps. 
To solve these relations, define the generating functions Pi(s;m) and Qi(s;m): 

Pi(s;m) =^p 1 (iV;m)s JV ; Q 1 (s; m) = £ qi (N; m)s N . (A3) 

2 

Eq. IA2|) then translates to 

P 1 (s;m)=8 m , + Q 1 (s;m)P 1 (s;Q). (A4) 

It is straightforward to find Pi(N; m) and Pi(s; m). The probability P i(N; m) is the standard random- walk proba- 
bility, given by 

1 r 27r 

P i(N;m) = — d0[cos9] N exp[im0]. (A5) 



From this we find the generating function (for m > 0) 



1 r 27r pimV 1 

Pl(s;m) = _L/ d6 e = _ (i _ s 2 ) i /2]r(1 _ 82) -i/ a (A6) 

Zir J 1 — s cos a s 
and the gambler's ruin (m = 0) generating function is 

P 1 (s;m = 0) = [l-s 2 ]- 1/2 . (A7) 
This yields for the gambler's ruin problem the well-known result 

Qi(s;m = 0) = 1- Pi(s;m = = 1 - [1 - s 2 ] 1/2 . (A8) 
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One learns from this that the probabilities for first return after K= 2, 4, 6, 8,. . . steps are 1/2, 1/8, 1/16, 5/128,. . . 
independent of the total number of steps in the random walk (if this number is larger than K) , and that the probability 
of ever returning is unity. This follows from Eqs. (| All IA3|I . which shows that Pi(l) diverges in d = 1 (and also in 
d = 2) but not in higher dimensions. 
For non-zero m we have 

Pi(s;m) =Qi(s;M)Pi(s;m = 0). (A9) 
It then follows from Eqs. (|X6l |A"71 lA"9)) that 

Qi(s;m) = {^[l-(l-s 2 ) 1/2 ]} m (A10) 

By expanding in s, one sees that the probabilities vanish for N < m, as expected, and by setting s = 1 in Qi(s;m) 
one sees that the probability of ever reaching m is unity. The probability qi(N;m) peaks for N ~ m 2 , when m is 
comparable to the vortex radius of gyration. 

Now consider calculating the probability of first return anywhere on an infinite plane in d = 3. The probability 
p 3 (N;rh) of going from the origin to lattice point m in N steps on an infinite lattice is given in Eq. i|14l and the 
corresponding generating function is 



P 3 (s;m) 

j 



i 3 r 2w s 

^[TT / d9,}{l- -[cos0i+cos0 2 + cos0 3 ]}~ 1 exp[i0-m]; P 3 (0;m) = 6 r% ff . (All) 

Take m 3 = in Eq. (|Allfl and sum from — oo to oo over mx t 2- The result is, not unexpectedly, a simple variant on 
the d = 1 gambler's ruin problem, and yields 

°° i fir fio 4 2 

«(.)- £ Pa(»; m,, mg , ms = 0) = - _ y _ - [1 - y + T ]-". (A12) 

1711,7712 = — OO u JO 

Using a standard expansion and Eq. I|A8|1 we find the probabilities qs(K): 



q3 (K) = -3- K ^C K 1/2 (^=) (A13) 



— 1/2 

where the C^- are Gegenbauer polynomials. 



2. Probability that a vortex pierces two separated Wilson loops [Sec. Iv| 

The next problem comes up in the compound Wilson loop estimates of Sec. where a vortex can pierce two 
separated Wilson loops. In d = 3, consider the plane surface of dimensions Li,L% centered on and perpendicular to 
the z-axis at a distance M along this axis from the origin. We ask for the probabilities Ps(N; Li,L 2 , M) to end up on 
this surface after N steps; the probability q^(N; Li, L2, M) of reaching this surface for the first time after N steps; 
and the probability of ever reaching it. Also to be calculated are the corresponding generating functions Ps(s; Li, M) 
and Qs(s; Li, M). We assume that all lengths Li,M are large in lattice units, so that the number of steps is also 
large. 

The first probability is a sum over the surface of the probability given in Eq. I|14fl : 

7Tli=Li/2 

p 3 (N;Li,M)= P3(N;m)\ m3=M . (AH) 

m<=-ii/2 

The sum is easily done to yield 

- ^&r* ] ^* + -* +m ^^* *%M n *°%M m - <ai5) 

In the large- TV limit N 3> L\, one finds an area factor emerging: 

A -A- r 27T 

p 3 (N;L u M)~ CT [TT / de j }[cose 1 +cose 2 + cose 3 } N exp[iMe 3 } (A16) 
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where A = L\Li is the area of the surface. This happens because the integrand is only appreciable when 9i < N^ 1 / 2 . 
The generating function for this probability, if needed, is constructed as in Eq. (|A1 1|) . including summing over to 1; m 2 . 
The probability q%{N; L\, L2, M) is determined by an analog of Eq. I)A2|1 . including a sum over surface variables: 

P3(JV;L i ,M)=5^53«fe(J;mi,M)p 3 (JV- J;m,,0) (A17) 

J rat 

where the sum over mx,m2 is delimited as in Eq. (|A14|) . To see what this sum means, write the ^-probability in 
Fourier form 

93 ( J\ mi, M) = J d 3 0q 3 (J; dy s ' n ~ l (Ai8) 

with ma = M. We find, using Eq. (fTljl . 

MN^M> . E / ^"Win^^ / ^[(E«-)/3]-. (Al 9 , 

We anticipate, and can confirm later, that in the limit N — > 00 the 9j are of order J -1 ^ 2 and the a-,- are of order 
(N — J)" 1 / 2 . It turns out that N — J is large compared to J (which is of order M 2 ) and so we can drop the ai.2 
in the argument of the sine functions in Eq. IjAlQJl . The a integral then factors out, and is given by the probability 
pz{N — J; 0) of returning to a given point after N — J steps. 

The problem is now solved, in principle, by using generating functions, as in Eq. (|A9|) . Or one can study the sum 
of Eq. (|A19(1 directly, and find by a scaling argument that the maximum value of the g-probability on the right-hand 
side behaves like A/ J ~ A/M 2 . Note that when lattice lengths Li, M are converted to physical lengths by multiplying 
by A this probability remains finite. 
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